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We derive an analytic expression for the zero temperature Fourier transform of the density-density 
correlation function of a multicomponent Luttinger liquid with different velocities. By employing 
Schwinger identity and a generalized Feynman identity exact integral expressions are derived, and 
approximate analytical forms are given for frequencies close to each component singularity. We find 
power-like singularities and compute the corresponding exponents. Numerical results are shown for 
O j N = 3 components and implications for experiments on cold atoms are discussed. 
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I. INTRODUCTION 

It is well known that as a result of spin-charge separation, interacting one-dimensional spin- 1/2 fermions with 
repulsive interaction at incommensurate filling form a two-component Luttinger liquid (LL). Multicomponent LL with 
more than two components can be obtained in fermionic systems with repulsive interaction^—. Possible realization 
of such systems are provided by multichannel quantum wires^, carbon nanotubes^^ and biased bilayer graphene 7 -. 
Interaction of acoustic phonons with spin-1/2 fermions in one dimension can also give rise to a multicomponent LLj^i 9 - 
In Mott insulating materials, a multicomponent LL can be formed in spin-orbital chains^— and spin tubesi 9 . under 
the effect of an applied magnetic fiel d 20 ' 21 and in spin-1 chains with biquadratic interaction s 22 ' 23 . More recently, 
atom trapping technology has permitted the realization of Bose-Fermi mixtures^! - — as well as degenerate gases 
with internal degrees of freedom . 28 ' 29 In the latter case, it has been suggested theoretically that these systems could 
realize SU(N) spin systems 3 ^— in low dimensions. In parallel, techniques for trapping atoms in one dimension 
have been developed 3 ^—. In Bose-Fermi mixtures trapped in one dimension, a multicomponent LL behavior is 
cxpccted*2r— . Similarly, multicomponent systems with repulsive interactions arc also expected to exhibit in one 
dimension a multicomponent LL behavio r 57 ' 58 The real space correlation function of the multicomponent Luttinger 
liquid can be readily obtaine d 8 ' 50 ' 59 . However, the majority of experimental observables are actually Fourier transform 
(FT) of these correlation functions and thus the FT of the multicomponent LL has to be obtained. Because of the 
branch cut structure of the correlation functions in a LL, this is a non-trivial task. In the single component case, 
the calculation can be done in closed former—. For two component systems, only the exponents of the power law 
singularities could be predicted^ 3 - - — Recently a closed form of the 2kp component of the density-density response 
function in a two-component LL (i.e. the spin-1/2 case with different charge and spin vclocites) at zero temperature 
was obtained^ 6 - in terms of Appcll hypcrgcomctric functions! 67 ' 68 Such an expression permits the description of the 
crossovers between the different power law singularities of the response functions. In the present manuscript, we derive 
an exact expression for the FT of the density-density correlation function in the general case of a multicomponent 
LL with different velocities for the modes. We show that in this general case, the FT of the Matsubara correlation 
functions are expressed in terms of the Srivastava-Daoust generalized hypergeometric functions. We give give the full 
analytical continuation of the correlation functions to real frequencies, recovering the leading power-law singularities 
and describing the various crossovers between them. 

The paper is so organized. In Sec. [TT] we give the Hamiltonian of a general multicomponent system and derive 
the general expression for the Matsubara correlation functions at zero temperature. In the subsections III Al and III Bl 
we derive two exact integral representations of the Fourier transforms of the Matsubara correlation functions by 
means of a Schwinger identity and a Feynman identity, respectively. The integral representation obtained from the 
Schwinger identity is used to predict the exponents of the power-law singularities. In Sec lII CI starting from the 
integral representation derived from the Feynman identity, we obtain the analytic continuation of the Matsubara 
correlator in various cases and give the asymptotic expression close to the singular points. Finally, we give some 
conclusions in Sec lIIII 
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II. MODEL 

In the present paper, we wish to consider the case of a general multicomponent system. The continuum Hamiltonian 



H = J^J dx 

^ / dxdx l V ah (x - x , )p a (x)p a {x'), 



(1) 



Ka<b<N • 



where ip a annihilates a particle of type a, p a {x) = tpt^a is the corresponding particle density, and M a the corresponding 
mass. The interaction of particles of type a with particles of type b is V a b- The particles may be either bosons or 
fermions. The Hamiltonian (jTJ) can be bosonized^ and its expression reads: 



/dr 
— [(7rn a )M a6 (7rn 6 ) + {d x <f> a )N ab {d m <i> b ) 
7T 



(2) 



where {4> a {x),Hb(x')] = i5 a b5{x — .t'), and the matrices M and N can be obtained respectively from the variation 
of the ground state energy with change of boundary conditions and particle numbers^ The Hamiltonian ([2]) can 
be diagonalized, and the equal time correlations can then be obtained^ We wish to calculate the Matsubara time 
dependent correlation functions. In the case of density correlations, since the density is expressed as: 



p a (x) = d x cf> a + A m cos[2m(0 a (x) - np°x)], 



(3) 



m— — oo 



We will need to calculate correlation functions of the form (T r e lm< ^ a ^' T 'e ml< Mo,o)^ wnerc y r [ s the Matsubara time 
ordered operator. If we have bosonic particles, the bosonized form of the annihilation operator being: 



Tpa(x) 



„i6 a {x) 



i2m[(f> a (x) — irp^x] 



(4) 



where V6* a = 7rn a , the leading term in the single-particle Green's function is proportional to (y r e l6 'aO,''') e - J M^ T )^ 
Because of the duality transformation^ M ■n- N and 8„ f) (t a , it is sufficient to calculate the correlation functions 

of the form (r T e'^" a "^' I ' T 'e" 1 ^»°"^' 0,0 ') and use the duality transformation to obtain the correlation functions 

(T Te i£„«nM*,T) e -i£„«n0„(o,o)^ Wc have . 



G nm {x,r)=^ / ^ e -l«l«(l-e^™ T ))[(^ + (MJV)g 2 )- 1 M] 7 



^ ^ ^ n & rn G n m ( X , t) 



2tt 



(5) 



(6) 



Let us introduce the projection operator P\ that projects on the eigenspace of eigenvalue u\ of the matrix MN to 
rewrite the matrix (a; 2 + (MiV)g 2 ) -1 as: 



(u;l + (MN)q 2 )-'=Y^ 



(7) 



In the limit of /3 — > oo, we have: 



+ (MNtfrHl e^™))e-l^ = £ ^(1 



Jqx-U\\qr\\ e -\q\a 



(8) 
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So that: 



and thus at zero temperature: 



g(*,t)4e^W * 2+m ; i+q)2 ), o) 



4 x ux 



2 \ V\ 

(rre iE„a^(^) e -iE»^„(0,0) )= -Q/ « ) , (10) 

- L - L \ + (ua|t| + a) 2 ' 



A 



with: 



»/a = > a a a b . (11) 

^ "A 

In the case of the correlation function, one obtains a formula analogous to (fit)]) , with the exponents r\ n replaced by 
fj n where: 

_ {Q\N) ab 

i "A 

a,o 



Q\ being the projector on the eigenstatc of NM having the eigenvalue u\. Wc note that for r = 0, using Eqs. (fTU|) 
(fTTj) . we recover the expression of the exponents of equal-time correlations frorn^. These formulas are particularly 
useful when implementing numerical methods like DMRG and exact diagonalization as one could use the following 
results to predict the response function from ground state energy computations. Knowing the Matsubara time ordered 
Green's function, Eq. (|10p . we wish to obtain the corresponding Matsubara response function: 

X M (q, iu n ) = J dxdTe l{qx -^ T \T T e l ^^ a ^ x ^e- l ^ aM0 ' 0) ), (13) 

and the retardated response function x(<7, ui) = \ M ili^n — > w + i0). In the following, we use two complementary 
approaches, the first one based on the Schwinger identity^ in Sec. Ill Al that will allow us to predict the singularities 
of the retardated response function, and the second one based on the Fcynman identity^ that will allow us to make 
a connection with the results for the two-component case^. We will assume that we can take the limit a — >• in the 
integrals, i. e. r\ = < 1. 

Let us remind that for the density-density correlation function its FT is related to the scattering cross section a of 
light at a frequency oj and angle Q incident on a sample, by the relation 

<x S(q,ui) = Sign(w)3mx(g,w) (14) 



where S(q, uS) is the dynamic structure factor. This quantity is accessible e.g. by means of inelastic neutron/light 
scattering when spin/density fluctuations are induced in the system and their subsequent relaxation is measured. 

A. Schwinger identity 

The Schwinger identity is2^: 

V Aa _ le _A(^+(„ T )=) (15) 



{x 2 + H 2 )« J Q T{a) 
In the multicomponent case, it allows us to rewrite the expression (|10j) as 



n ( ^ ,(,, m + ^ ) /n r(2\ exp 



-x 2 



^x 2 + {u„ \t\ + a) 2 J J n r (^«) 
The Fourier transformation in (|13[) reduces to a Gaussian integral, and we find 



Em- t2 £ A -< 



(16) 



X (q, ,:„,„) = ng^vn / TT " " e / L (17) 
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With the change of variables: 
we rewrite (fT7|) as: 



Xj = = 2, . . . ,N), 



(18) 



r(ryi) 



N-l j ..fi + i- 1 



, . x^N — l n 2 l v^^V— 1 2 



r ( r /j+i) //-. . ^jv-i \ 7 2 , v^jv-i a" 



(19) 



With the change of variable Ai = /j 1 , the integration over Ai can be done in closed form. We find, provided that 
EVn < 1: 



X(g,iw„) = ?ra 



2El- 



n.rfe) 



iV-1 



,»)j+l-l 



, ; -l 1 2 , v^jv -i 5 — 



(20) 



We can perform the analytic continuation on (|17p by substituting iw„ — > lo + iO. When oj ~ u\q, we have to consider 
the integral: 



/'nW- 1 1 



2 2 2 i ^22 2\ 

j=i 



(21) 



With the change of variables /jj = (ufq — 1 we finally find that for |cj| — > uiq, 

x(q,u)~\(u iq ) 2 -u 2 \ 2 ^^-^-\ (22) 

provided 2(^? / , „) — 771 — 1 < 0. In the general case, we expect to find x(<Z: w — u jq) \( u jq) 2 ~ uJ 2 \ 2 ^ n ''"^''j -1 
when 2(Y] Vn) — f}j — 1 < 0. We note that 77 = is the exponent of the equal time correlation functions.— The 

exponents of the singularities satisfy a "sum rule" : 



N 



J2(2v-Vn-l) = (2N-l)r)-N 



(23) 



B. Feynman identity 

Using the identity fromSi (Appendix B), we can rewrite the correlation function (|10p as a multiple integral: 



N TVY^ 2 \ r N N 

f} i (x 2 +u 2 n r 2 )^ n?=ir(jn) J ji 3 U 



N 



-(£f=i"i) 



3=1 



3 u, 2 r 2 



(24) 



The Matsubara response function (fT3"|) is then found from the intcgral(11.4.16) of Ref. UM 



J dxdre 1 ^-^ {x 2 + Y^ w jU 2 j t 2 ^ 

22(1-7)^(1-7?) / 2 , 



JV 2 



(25) 
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where 1/4 < r) < 1 as: 

p/i s /• N N 92(1—17) 2tj I n \ 




where we have defined 77 = J2j Vj- For the case of 77 = 1/2, the integral (f2"6")l can be expressed as a Lauricella function 
Fo, which actually reduces to a simple product. The result is simply: 

N 

X (q, iu> n ) = 2na [](u; 2 + utf)-* . (27) 
3=1 

This result generalizes the one obtained for the two-component case in Ref. [66l 
When r\ > 1/2, we can use again the Feynman identity^ to write: 

^ - ix, - - J ' J8S " (1 - a) '" /2 f V + g *») ' ' (28) 



(#+E J «H " (E. 

and find: 



(M) . . 7 r2 2 ( 1 -") a 2 r(i/2) ^ 2( 1} f 1 s -^i-^-3/2 



p(iV-l) f 1 "AT - ~ U N-1 , on , 

F D l-; Vu ...,r lN . 1 ; m 2 — Z (29) 

where the label (M) stands for the Matsubara correlation function and i 7 ^ is a Lauricella hypergeometric func- 
tion of TV — 1 variables. The integral (|2"9")l can be expressed in closed form^ using Srivastava-Daoust generalized 
hypergeometric series^ of N variables as: 



(1/2; 1,..., 1,1) : 771, . . . ,7/jv-i, 1 - 7/ «? <-i u 2 

(77;!,. ..,1,0) : -, -, 1 u f v '---' 1 M ^ ' u 2 N q 2 



(30) 



where F(?y) = x2 - p^J^ 1 — ■ When N — 2, the Srivastava-Daoust hypergeometric series reduces to an Appcll -F2 

hypergeometric function of 2 variables. Using the identity (16.16.3) from [68|, this function is seen to reduce to the 
Fi Appell function expression given in Ref. [661 ] . Since we are not aware of any study of the analytic continuation 
of the Srivastava-Daoust hypergeometric series outside their circle of convergence, we will not pursue with Eq. (|30[) . 
Instead, we will consider directly the analytic continuation of the integral (|29|) . Introducing the new variable: 



the integral (|29|) is rewritten: 



; = 4±44 s , (3i) 

sw n + U N<1 



r(?7 - l/2)r(7?)MAr \ 11% J Jo 

(JV-i)/l. . _u 2 N -u\ ( U 2 n t \ _ U% - U%_ 1 U%t 



Fd [ ^i 7 ?!) • • • i^iv-ii 7 ?; -g-_zi ( 1 _ )..,,; g 2 -M l- I , (32) 

V 2 <% V ^n+ u WJ u n V ul + u z N q 2 J) 

In the case N = 2, the Lauricella function in the integral (|3"2"|) reduced to a 2F1 Gauss hypergeometric function. Doing 
the integral gives a F\ Appell hypergeometric functional. With N = 3, the Lauricella function reduces to a Fi Appell 
function. 
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In the case of 77 < 1/2, the integral (|29|) is divergent. In order to obtain a convergent integrals, we write in Eq. (|26[) 

n-i 



w + 1 Ej "./"'.. 



TV 



E, u>j 



7,-1/2 







1-7, 



and apply the Feynman identity to each term in the sum. We thus find: 

r(3/2) 



Ej 



X (g,io; n )=7ra 2 ^ 2 ''- 2 2 2 ( 1 -'') 



r(jj + i/2)r(77 + i) 



/' 

.7 



ds- 



3/2" 



JV 



2 F (JV-1) 



+6 ji } 1 <j<tf- 1 ;T) + 1; 



£=1 



S-^T + W 



Each term in the sum is then expressible with Srivastava-Daoust hypcrgcometric functions 



' • ■ • ' ,.,2 



(33) 



(34) 



C. Analytic continuation of the Matsubara correlator 



To obtain the retardated response function, we have to find the analytic continuation iio„ — > to + ie of (|32|> . We will 
first discuss the special case of 77 = 1/2, where the continuation is straightforward, leading to a simple picture of the 
behavior of the retardated response function. Then, we will turn to the more complicated case of 77 > 1/2, for which 
the calculations are more involved. We will see however that the simple picture of the case 77 = 1/2 is preserved. 



1. The case of 77 = 1/2 

In the case of 77 = 1/2, the analytic continuation is easily obtained from Eq. ([27)) . We have for ujq < ui < Uj + iq: 



N 



X(q,u + i0) = 2iTae m ^ 3 '=^' u JJ \lu 2 - ( Uj q) 2 \-^ 

i=i 



(35) 



So that: 



3mx{q,u — > Ujq + 0) ~ 27rasin 



1=1 



and: 



3mx(q, to — > Uj + iq — 0) ~ 2ita sin 



L 1=1 



(«i5) a r % ni(«i- u iVi-'"> 



¥j+i 



(36) 



(37) 



showing that the spectral function has singularities with exponent —rjj everytimc to ~ w^g. This result is is agreement 
with the result of Sec. Ill Al for r\ = 1/2. The spectral function has a threshold for lu < u\q as there are no excitations 



(38) 



which implies that the peaks of the spectral functions are asymmetric around ui = Ujq. The imaginary part of the 
response function is represented on Fig. [T] for the case 77 = 1/2. 



Moreover, we note that 


for 


7>1: 


3mx(q,u -» u 3 q + 0) 


sin 


t E;=i m 


3mx(<?,w -> u i1 - °) 


sin 


^EiU m 
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FIG. 1: The imaginary part of the response function x(q,iS) for the case rji = 0.25, 772 = 0.15 and 773 = 0.1. The velocities 
are ui = 2ui and U3 = 3ui. The unit of frequency w is uiq. The unit of 3tnx(?, u) is 2na/(uiq). Power-law divergences are 
obtained for uj = «i,a,39. 



2. The case of rj > 1/2 



In the general case, we need to find the analytic continuation of Eq. (|32p for 77 > 1/2 or (|34p for 77 < 1/2. Since the 
method is similar in the two cases, we will concentrate on the case of 77 > 1/2. 
Formally, under the analytic continuation, the variables in Eq. (|32j) become: 



'JV 



N 



1 - 



+ u 2 N q 2 



1 




(w + ie) 2 



(it/vo 1 ) 2 - (u + it) 2 



(39) 



The Lauricella function Fd bas cuts every time the real part of one of the variables is larger than one. In the case of 
lu > ujvQ, the real part of all the variables, acconding to (|39|) will remain less than one, and the analytic continuation 
is straightforward. We can then derive an equivalent for the Lauricella function in the limit of ui — > u^q + 0. We find: 



F 



(JV-l) 



D 



1 

2 57?1 " 



N 



'N 



1 - 



9 2 9 

uj- - u z N q z 



N 



'JV-l 



)*•*)! 



1 - 



UJ Z t 



- u 2 N q 2 



r(77)r(7 ?A r + 1/2 - 77) /w 2 - (u/vg) 



r(i/2)r(r ?JV ) 



2 \ v-vn JV— 1 



0J 2 t 



JV 



->i., 



yielding: 



JV-l 



= 7 r2 1 -^ 2 r(77 jv+1 - 277)sin(^) (oj 2 - (^g) 2 ) 2 "-^- 1 , i 



(40) 



(41) 



In Eq. (|4"Tj) . the exponent predicted in Sec. Ill Al is recovered. For u 
estimation of the Lauricella function in the form: 



(*)' - 

uj^q — 0, it is also possible to find an asymptotic 



F 



(Jv-i) 



D 



1 



N 



'N 



1 



9 9 9 

u - u n1 



2 2 

U N ~ U N-1 



1 ■ • ■ 1 7 



N 



0J 2 t 



9 2 9 

oj z — u N q z 



r(77)r(77 A r + 1/2 - 77) f(u N q) 2 - 



T(l/2)T( m ) 



2\ V-VN 



ujH 



JV-l 



i7r(r)-r) N )sign(w) 



(42) 



giving: 



JV-l 



3mx(g,u; -> tj/vq - 0) 



7r2 1 -"a 2 r(77 / y +1 - 277) sin[7r(77 - 77^)] (( Ujvg ) 2 - q^-w-i 

T(r) N )u N . ~. f..* W- 1 11 I 



3=1 



-Vj 



(43) 
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i. c. the same power law divergence as in Eq. (|41[) but with a different prefactor. The ratio of the two expressions is: 
sin[7r(?7 — ?7iv)]/ sin(7T77), as previously noted in the special case of 77 = 1/2. For e — > 0, the imaginary part of (j3T))) is 

positive for w < uj^q and the real part of (l39| is equal to one for < = U = \^ Nq /, ul \ 2 Z\ ■ We have t\ < t 2 < ■ ■ ■ < ijv-i. 

In particular, one finds that when uj < Uiq, t\ > 1, so that no analytic continuation of the Lauricella function under 
the integral sign in ()32|) is needed. The response function remains purely real in that case, and the spectral function 
vanishes. In the case Ujq < uj < Uj+iq, we find that t\ < . . . < tj < 1 < tj+i < . . . < tjv-i- As a result, the integral 
has to be split into a sum of integrals over the intervals [0, t\], [ti, tj+i] with 1 < I < j — 1 and [tj, 1]. For each interval, 
the analytic continuation of the Lauricella function Eq. (|B3[) must be used in order to express the full integral. To 
give a concrete example of the procedure, we will at first focus on the case N = 3. We need to consider the following 
integral: 



l(q,cj) = [ dtt-^l-tr- 3 / 2 F 1 (l; m ,r ]2 - 1 r ] ;^ 

JO \ z u 3 



1 - 



u\q 2 



1 - 



u\q 2 



(44) 



For lo < U\q, we have: 



Jo V 2 u 3 V W39 W / u 3 V W3S w // 

Let us now consider the case of Mig < w < First, we split the integral with the rule given above and consider 
the intervals i g [0, — t ] and i G [ — "5 , 1]. In the first interval, the analytic continuation is straightforward. 

— 1~ 1 — T — 1 
u l "l 

In the second integral, we must use (|A4[) . 

To calculate the integrals, it is convenient to perform the following change of variables in the first and second 
interval respectively: 

t = -3-3--! 



Ul 



U3 ^ I " 3 \ 



S2 



In the integration over si, we don't need an analytic continuation of the Appell function. In the S2 integration, one 
of the two variables is larger than 1 and we find the analytic continuation using Eq. (|A4[) . 

Since we are interested in calculating the imaginary part of (|44|) . only the term proportional to e I7rj ' 1 gives a 
contribution and we can perform explicitly the integral when \oj\ — > U\q. The final results for the imaginary part of 
the integral (|4"4")l yields: 



3ml(w, q) = 



r(l/2)r(l/2 + ry- )?1 ) /^jn^-i 




dss „- ( , 1 -l/2( 1 _ s )r,-3/2 / j j W 2_ u 2 2 2 ^2 ^2 _ ^2 2 

*i (i - 2^2;r? + - - m; , , 1 J J.. 2 _ 2 1 2 * )(46) 



When w — > Mig, the integral is behaving as (to 2 — g 2 ) 2j?_1_,?1 , in agreement with (|2"2"]l . For w — > we need to 
consider the behavior of the Appell function as one of its arguments is going to unity, while the other is negative. 
Using the results from App. [C] we find that when 771 < 1/2, we can use Eq. (|C8[) to approximate Eq. (|4"o]) as: 



1 /„.2 „,2 C„. „\2 ,.2 \ V-m-l/2 



7mI{ U , q ) oc / ^f V 8 !^ 2 "" 

Jo \ u 2 ~ M l (""3?) - W 

I - ?J 2 (7J29) 2 - W 



9 

2 \ 2))-?) 2 -l 



uf (U39) 2 - ^ 2 



(47) 
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Again, this result is in agreement with the power law divergence expected from (|22|) . 

We now turn to the case u 2 q < oj < u 3 q. In that case, we have to split the integral in Eq. flH| into three 

integrations. The first one, on [0, ] does not contribute to the imaginary part. The second one, on the interval 

[ (ua/tt^j 2 -! ' ("a/tttj 2 -! ]' recm i res an analytic continuation of the F\ function using (|A4[) and gives a contribution to 
the imaginary part: 

r(i/2)r( m )r(i/2 + 7j-m) WJ U§-^17 V«2-«?/ U 2 -MV V 



^3 - U 2 J V M l M 3 - «2 / V U 3 - U 2^ - [V-lVr 

( 1 1 ul - u? \ 

xFi 1-771; + ^ -J71; — 2 2 s ' s ( 48 ) 

The last integration, on [ ("'^"^ _ 1 , 1], requires an analytic continuation of F\ using (|A7|) . and contributes two 
terms. The first one is: 

7Tr(7?)r(l - r ]2 )\w\(u 2 3 ) 2 ^^ 1 (u 2 2 )- r !(u 2 1 )^ + 1 ^ (u 2 - M 2)r,+ l/2-2, U -r, 2 ^2 _ (^^2)^-1/2 

r(i/2)r(?7i)r(^ - i/2)r(i -m- m) («§ -ui) 1 -'' 1 -'' 2 ^ - u 2)-,-m-i (( M3g ) 2 - ^y> 

x *(! - s)-3/ 2 ( ! + n W-{u 2q f \ " / ^-(^ \ ^ ( u 2 -uW- { u 2 q) 2 
Jo [ j V uKu.q) 2 -^ 8 ) [ L+ (u a q) 2 -u; 2S ) + u\ - u\ (u 3 q) 2 - UJ 2 



13 u\ — u\ 1 

Fx | 1 - Vl ; -, - - q; 2 - m - y 2 ; -g— g, \ , (49) 

"'2 ~~ U l ( u '^ ( i) 2 " uj2 



and for 77 — 7/2 — 1/2 < behaves as {uj 2 — (i£2<?) 2 ) 2?? m 1 in agreement with ([22]) . The second one is: 

rfa)r(l- 772) sin[7r(»ft+ 772)] /^- U 2\"Y^ 2 u\-u\ y { uju 2 - (u 2 q) 2 Y>~ 1/2 fu 2 ~ (u 2 q) 2 Y>~ V2 - 1/2 



u 2 (u 3 q) 2 — lj 2 ) \cu 2 u 2 — u\ ) \(u 3 q) 2 — u 

/ 1 1 a; 2 — (u 2 q) 2 ul — u 2 lu 2 — (u 2 q) 2 \ , , 

x^ 1 - ry 2 ; -, r? i; r? + - - W - ? — \ h ^^2* 50 ) 

V 2 2 ("3?) - u^- u{ (u 3 q) 2 - lu 2 J 



The latter term contributes a divergence {lu 2 — (u 2 q) 2 ) 2v ^ r ' 2 ~ 1 as w -> W2<Z + in agreement with (|22|) . 

To summarize, the qualitative behavior of the spectral function is the same as in the special case of rj = 1/2. The 
spectral function has a threshold at uj — u\q, and has power law singularities for uj = Ujq with an exponent given by 
Eq. (EH). 

The case of a general TV is treated in the Appendix [D] It is found that the leading singularity for lj — > Ujq is again 
a power law divergence, with exponent given by (|22j) . 



III. CONCLUSION 



We derived a analytical expressions for the zero temperature Fourier transform of the density-density correlation 
function and the bosonic Green's function of a multicomponent Luttinger liquid with different velocities. By using 
both a Schwinger identity and a generalized Feynman identity, we derived exact integral representations while an 
approximate analytical form was given for frequencies close to the characteristic frequencies of the different collective 
modes of the system. We derived in detail the analytic continuation for generic and discussed, as an example, 
the case TV = 3. Power-law singularities arc found every time the frequency is equal to the characteristic frequency 
of a collective mode (u}j(q) ~ u^q), with the same exponent 2r/ — rjj — 1, but a different weight when approaching 
the singularity from the left or from the right. The power-law singularity replaces the expected delta function for 
nonintcracting particles. Moreover if the characteristic exponent at the singularity becomes negative, as in the case of 
systems with attractive interaction when considering density-density correlations or in systems with repulsion when 
considering bosonic Green's functions, a cusp is expected to replace the power-law divergence. All these results are 
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FIG. 2: Spectral function close to the singularity for N=2 and the values of the parameters: r/i — 0.3,772 = 0.1,77 = 0.6, ui = 
0.25, u 2 = 0.75. 



valid in the ground state. For nonzero temperature T, the power-law divergence at uj = Ujq is replaced by a maximum 
diverging as T 2v ~ Vi ~ l as T — > 0. In the vicinity of the maximum, from simple scaling, we expect Imx(g,w) ~ 
(fc_B? 1 ) 2l '~'' i ~ 1 /tcmp.[(w — Ujq)/kBT], with /temp. (a; 3> 1) ~ aj 2 * 7- ^ -1 . If we now turn to a zero temperature system of 
finite length L, we expect the power- law divergence at u = u^q to be replaced by a maximum diverging as (1/L) r, ~ ri ^ , 
and in the vicinity of the maximum Imx(q,oj) ~ (1/ L) 271-71 ^ 1 fi en \L(uj /uj — q)] with /i on . ~ (x ^> 1) ~ x 2r >~ r >j- 1 . If 
the calculation of both functions /i on . and /temp, remains an open problem, the scaling arguments suggest that the 
power law behavior of Imx(g l , u>) is observable for finite temperature and finite size provided \lu — ujq\ 3> fc^T, Uj/L. 
Such behavior could be probed by Bragg or time of flight spectroscopy in the case of atomic gases of mixed species or by 
inelastic neutron scattering technique in one quantum magnets with orbital and spin modes. An obvious extension of 
the results of our manuscript is the calculation of fermion (or anyon) spectral functions in multicomponent Luttinger 
liquid. In the case of a two-component liquid, the fermion spectral functions are expressible in terms of Appell 
hypergeometric functions^. For the case of three or more components, the results of our manuscript hint that 
the fermion spectral functions should be expressible as Srivastava-Daoust hypergeometric functions or a suitable 
generalization. 
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Appendix A: Analytic continuation of the Appell Fi function 



In the calculation of the response function, an analytic continuation of the Lauricella hypergeometric function 
^ is necessary. In the present appendix, we present the analytic continuation of the Appell F\ hypergeometric 
function that corresponds to the particular case of N = 3. The analysis of that case is the stepping stone for the case 
of general N. 

In order to find the analytic continuation, we start from the integral representation of the Appell hypergeometric 
function£!££ 

r(c) f 1 t a_1 (l - rj) c_a_1 
F 1 (a;fe 1 ,& 2 ;c;, 1 ,, 2 ) = r(a)r(c _ Q) y o dt ^—^—^ (Al) 

when \zi\, \zs\ < 1, expanding in series (|Aip and integrating w.r. t. t gives back the series expansion. We wish to 
use (|A1[) to express lim £li£2 _,.o + F\{a\ b\, 62; c; Xi + ie\,X2 + K2) for x\ and x-i real for the cases x-i < 1 < x\ and 

1 < X2 < Xi. 



1. case X2 < 1 < xi 



In order to calculate the integral (|A1[) , we need to take into account the sole branch cut of (1 — t(xi + ie\)) bl for 
t > \jx\. Because of this cut, we have for t > xi, (1 — t{x\ + «ei)) bl = e~ mhl {tx\ — l) bl . Therefore, we split the t 
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integral in (|A1[) into two integrations on [0, 1/xi] and [1/xi, 1]. After a change of variable t = s/x\, we find for the 
[0, xi] integral : 

1/xi t a ~ 1 (i-t) c ~ a ~ 1 i r(a)r(i — £>i) 

d \l-t Xl )Hl-tx 2 )»> = ^ T(l + a- bl ) Fl{a; « + 1 - c, 62; a + 1 - 6i; x 2 /xi), (A2) 
a purely real expression. Then, for the [1/xi, 1] integration, we find: 

J 1/Xl (l-t(x 1 +z0+))" 1 (l-te 2 )^ xr b ^ 1 (x 1 -x 2 )^ r(l + c-o-6i) 

xFi ( 1 - 61; 1 - a; 62; c - a + 1 - 61; 1 - xi, a?a ^ a?1 ~ 1 M , (A3) 
V xi - x 2 / 

where we have used the linear change of variables t = 1/xi + s(l — 1/xi). Our result for the analytic continuation is 
then: 

Fi(a\ 6i,& 2 ;c;xi + i0 + ,x 2 ) = — ^-^77 , — , , ^F 1 (a;a + 1 - c,6 2 ;a + 1 - 6 X ; l/x x ;x 2 /xi) 

1 (c — a)l (1 + a — Ox J x" 

„ fei r( c )r(i-6i) (x! - i) c -°- bi / ^(xr-i) 

+e 1 : , , TT-^b^TT HT F i 1 - 61 1 - a 6 2 c - a + 1 - 61 1 - xi, 

r(a)r(l + c- a - 61) ^ b2 L ( Xl ~x 2 ) b2 V xi-x 2 



2. case 1 < xi < x\ 

In that case, we need to consider both the branch cut of (1 — t(x\ + i0+)) bl and of (1 — t(x 2 + i0+)) b2 . Thus, we 
are lead to split the integral (|A1[) into three integrations over [0, 1/xi], [1/xi, l/x 2 ] and [l/x 2 , 1]. The first of these 
integrals is still given by (| A2|) . The second integral is given by: 



1/3:2 t^ji-ty-*- 1 = ivbi r(i - b 1 )T(i - b 2 ) (xx - X2 ) 1 - b ^-(x 1 - ly-*- 1 
i /xi (i-t(x 1 + io+)) b i(i-te 2 ) b2 e r(2-6!-6 2 ) x c 1 - b *- 1 x 1 2 - bi 

XFi ( 1 - 61; 1 - a, 1 + a - c;2 - h - 6 2 ; 1 - ^, 7 ~ ^ 1 , (A5) 
\ x 2 x 2 (xi - 1)/ 

where we have used the change of variable t = 1/xi + s(l/x 2 — 1/xi) and taken into account the branch cut of 
(1 — t(xi + iQ+)) bl . For the third integral, we have: 

1 dt t^ji-ty— 1 = cl7r(bi+b2) r(i-b 2 )r(c-a) (x 2 ^ir^ 

1/X2 (1 - t{ Xl + M+)) b i (l-t(x 2 + i0+)) b * r(l + c-a-6 2 ) x^ 6l (xi -x 2 ) b i 

xFt ( 1 - 6 2 ; 1 - a, 61; 1 + c - a - b 2 ; 1 - x 2 , ^C 1 "^) ] ( A 6) 

V xi - x 2 y 

where we have taken both branch cuts into account, and we have used the change of variables t — l/x 2 + (1 — l/x 2 )s. 
The final result is: 

Fi(a, 61,62; c; xi + i0+,x 2 + i0 + ) = ~~^7T ~ r ; A-^ 1 ( a; 1 + a - c, 6 2 ; 1 + a - 61; — , — 

1 — djl ( X — (— GE — Ox) Xi V Xi Xi 

, F ! C J r (1 - W 1 ,- b2 ! ; ( ^ = 1)C ' a " f 1 - hi 1 - a, 1 + a - c;2 - 61 - 6 2 ; 1 - * J* 



r(a)r(c-a)r(2-6i-6 2 ) x c-b 2 ~i x i- bl ^ - . . *. ^'^(^ 

, h 1 r(c)r(i-6 2 ) (x 2 -i) c - a - b2 / , xi(i-x 2 ) 



Pi 1 — 69; 1 — a, 61 ; 1 + c — a — 69; 1 — x 2 , ■ 
r(a)r(l + c-a-6 2 )4- b i( Xl -x 2 ) b ! V a* - a* 



For the case AT = 3, the expressions (|A4[) and (|A7|) must be injected in the integral (|3"2"|) after analytic continuation 
to yield the response function. 
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Appendix B: Analytic continuation of the Lauricella Fn function 

The Lauricella Fd function has the integral representation: 

r(c) 



F D (a;bi,...,b n ;c; z 1 ,...,z n ) 



1 t a-l( 1 _ t }c-a- 

r(a)T(c-a)J *m=i(l-^) 



If we want to calculate for x n < x„_i < . . . < X2 < X\: 

lim F D (a; 61, ... , b n ; c; xi + ie±, ...,x n + ie r . 



(Bl) 



(B2) 



we have to consider the cuts of the functions (1 — txj + iQ^)~ bj . We are thus led to consider separately n cases, i. e. 
x n < . . . < Xj+i < 1 < Xj < . . . < Xi with j = 1, . . . , n — 1 and 1 < x n < . . . < x%. 

In the case of x n < . . . < Xj+i < 1 < Xj < . . . < x%, the integral (|B1[) has to be split into j + 1 integrations over 
the intervals [0, l/xi],[l/xi, I/X2],. . . ,[l/xj, 1]. In analogy to the case of the Appcll F\ function, each integral gives a 
contribution proportional to a Lauricella Fd function. The resulting expression is: 



F D N '(a;b 1 ,...,b N ;c;xi + iO, . . . ,x N + iO) = 



r(c) 



F(a)F(c-a) 
3-1 



r(o)r(i - h) f (n) 

T(l + a-h)xl D 



a: 1 + a - c, 62, ■ • ■ , b N ; 1 + a - b 1 : 



1 x 2 



E 



r(l - b m )T(l - bm+i) (x m - x m+ i) 



l-b m -6„ 



+1 (l^m) C 



' x L(2-6 m -6 m+ i) 



c-l-6 m+ i l-6 m 
I'm x m+l 



F b N) 1 _ b rn, bi,..., b m -i, 1 - a, 1 + a - c, 6. 



n 

1 - 



Xi Xi 
X/ x^ 



XjV 
Xl 

61 



7 



m+2, •■• , »JV, 



b N :2 - b 



1 _ £nL 



, 1 



1 



Xm+l 1 X^ 



3? m + 2 



r(l-6 3 -)r(c-a) (x, - 1) 

ra + c-a-b,-) 



c—a — bj 



«c-l 



n 



IE/ 



-61 



(A r ) 



1 - x 



1 - &j-;6i, . . . 1 - a,bj +1 , . . . , b N ; 1 + c - a - _ ^ 



1 -x 



1 — r • " J 

Xj-! X j + 1 



1 — X; 



1 — I, 



(B3) 



By reducing to N = 2 it can be checked that the results of Sec. O are recovered. 

In the case 1 < x\ < . . . < x„, we have to split the integral (|B1I) into n + 1 integrations over the intervals 
[0, 1/xi],. . . ,[l/xj, 1/xj+i],. . . ,[l/x„, 1]. Each integration contributes a term proportional to a Lauricella Fd function. 



Appendix C: Asymptotic expansion of the Appell F\ function 

We wish to obtain an asymptotic expansion of the Appcll function Fi(a\ b\, b 2 \ c; xi, X2) in the case of X2 < and 
x\ — > 1_. First, we need to obtain an expression of the Appell function in the form of a convergent series for all 
X2 < 0. We consider the series expansion for the Appell hypergeometric function F\: 

ri / t , \ ( a )n 1 +n 2 (bl)n 1 (b2)n 2 x" 1 X^ 2 , . 

F 1 (a;b 1 ,b 2 ;c- 1 x 1 ,x 2 ) = > j-t r — r , (CI) 



(c) ni +n 2 n\\ n 2 \ 

(ffl)ni(6i)m ^i 1 



where we have used the notation^: 



£ V 7; 1Ml ^.2*1 (a + m, 62; c + n i; x 2 ), (C2) 



r(n + a) 

(a) - = (C3) 
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Using the second line of (|Cip . we can define the function F± for all x 2 £ [1, +oo[. For x 2 < 0, we can use Eq. (15.3.4) 
from Ref. 0, to rewrite: 

Fi(a;b 1 ,b2;c;xi,x 2 ) = (1 - x 2 J 2 2^ p F i I &2> c - a; c + m; — — 

n \ c )ni n l- \ X 2 I 

= (i-x 2 )~ b2 ^2 *- a ' > " i ^ i ^" i ^ 2 ^" 2 *- c ~ a ^" 2 x " 1 1 ( X2 



(c) 



(1 — X2) h2 ^3 ( a,c - a; 61, 6 2 ; c\x\, 



m! n 2 ! V^2 - 1 

■i'2 



a?2 - 1/ ' 



6 2 (c- a)„ 2 (fe 2 )n 2 / £ 2 



(c)„ 2 n 2 ! \ x 2 - 1 



"2 



2-Fi(a, 61; c+m;x-i) 



(C4) 
(C5) 



where we have used the convergence of the Gauss hypergeometric series to obtain the last two lines, with this 
convergent series, we can analyze its behavior as X\ — > 1_. First, when c + to — a — b\ > 0, we have from Eq. (15.1.20) 
in Ref\T& 



r(c + m)r(c + to — a — 61) 

hm 2^1(0, &i;c + m;xi = — r-^. 

z->a- 1 (c + m — ajl (c + to — oij 



When c + m — a — bi < 0, using first Eq. (15.3.3) in Refill 

we find that as x — > 1_, 
-0-61+m r ( c + m ) T (a + h-c-m) 



2 F l (a,b 1 ;c + m;x 1 ) ~ (1 - x x ) 
Thus, when c — a — b\ < 0, we have that: 



r(a)r(6i 



F 1 (a; h; b 2 ;c; x u x 2 ) ~ (1 - ^^ ^^ C) (1 - Xi)—^ , 

r( a )r(6i) 



while for c — a — 61 > 0, 



lim Fi(a;br,b 2 ;c;xi 

Xl— >1_ 



, , . r(c)r(c+o-6i) / 

T(c- a)T(c-h) \ 



i,6 2 ;c- &i; 



•'2 



012 — 1 



Let us now return to the case of a general N. We have to consider the integral: 



i {N - 1} ( q ,u)= f'dtt-^i-tr-^F^ [ i{ Vj } 



jn<j<N-i\r]] 



N 



1 



(C6) 



(C7) 



(C8) 



(C9) 



(CIO) 



l<j<N-ls 



Appendix D: Analytic continuation of the Matsubara correlator for generic r\ and N components 

With the help of Eq. (|B3I) we can write for ti < £ < £/+i : 

t 



where: 



F, 



(N-l) I 1 r i I "at — «j 

1 ^{vj}i<j<N-i;m > 



1-1 



D 



1 



N 



"81^ _ 1 , 



= £>„,(*) + Vi(*)> (Dl) 



T( V )T(1- Vl ) 



r(7?-i/2)r(3/2-TO) 
/ 



1 1/2 



(u+ioy; 



F 



(JV-1) 



13 3 



V 



2 9 



2<J<JV-1 



(D2) 
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for m > 1: 



rfa)r(i - r? m )r(i - r? ; 



'm+lJ 



'm+l "m 



2 \ !-»/»> 



r(i/2)r(T? - i/2)r(2 - Vm - - Vm+1 ) y u %- u 2 m+1 

n-3/2 



l m+l "m 



N u m 



2 \ -Vm+l 



n 

fc^m, m+l 



XT >7* 



1-1 



1 + ^ 



1 - 



— ^iT ^ ( 1 _ 7 hn, {Vk}l<k<m-1, 9' 2 _ 7 '' {^}™+2<fc<A r -i; 2 - 7 7™ _ ^m+l! (D3) 



+1/ J Kfc<r 



( U m U m+l) 
-1 ( U JV — U m+l) 



« -«fc)( U AT - M m+l) 



and: 



r(l - 7?;)r(77 - 1/2) 

rfa - ^ + 1/2) 



n 



-<//, 



'Ei »7fe 



feA"JV "m+i; J m+2<fc<AT-l 



■O-rH-1/2 



"W 



-, 7,-1 



(D4) 



J) [ 1 - rji, {rjk}i<k<i-i, o) {%}i+i<fc<Jv-ii ~ 



JV 



1-1- 



"N I 



N 



u ~n1 



Kk<N-li 



So that for ujq < uj < Uj + \q: 
l- 1 rti+i 



l-l 



\m=0 



'3-1 



\m=0 



2 2 

ufj - uj 



(D5) 



l<j<iV-l, 



We can rearrange that sum into 
3-1 r i 



(D6) 



m=0' ,I "'+ 1 



1=1 



+ I dtt-^i-ty-^F^-^ ( ^,{rii}i<i<N-i;rr, 



9 2 

< - uj 



N 



"iVT 1 
— 73 1 . 



l<3<JV-Xy 



So we have to calculate 2j integrals. When in the last integral we do substitute the expression of (|D4|) . a factor 

[ui 2 — (ujq) 2 ^) 2r> r>3 1 appears with an integral which is regular in the limit oj — > UjQ, so the asymptotic behavior 
previously predicted is recovered. The first class of integrals in (|D6[) can be, instead, manipulated by using the 
definition (|D3j) and performing the following change of variables: 



/ {u N q) 2 



1 



/ (u N q) 2 



1 



(S- 1 ) V ( 



(D7) 



The integrals turn out to be regular for ui > ujq and for ui — > Ujq + they behave as (uj 2 — {ujq)' 



giving only 
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a subdominant contribution. 
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